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Abstract. The technical procedure of a search 
of bursts of very high energy gamma rays from 
evaporation of primordial black holes on air-shower 
array "Andyrchy" of Baksan Neutrino Observatory 
of Institute for Nuclear Research is described. The 
theoretical model used in the present work assumes 
that the chromosphere around the evaporating black 
hole does not form. For minimization of the cosmic 
ray background the method of multidimensional 
analysis of modelled as well as experimentally de- 
tected events is appUed. The new upper limit on 
the concentration of evaporating primordial black 
holes in the local region of Galaxy is obtained. The 
comparison of the results of different experiments is 
given. 
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I. Introduction 

Primordial black holes (PBHs) can be formed in 
the early Universe through the gravitational collapse of 
primeval cosmological density fluctuations. Therefore, 
the formation probability of PBHs and their observa- 
tional manifestations depend significantly on how the 
primeval density fluctuations emerged and developed. 
Theoretical predictions of the PBH formation probability 
depend strongly on the adopted theory of gravitation and 
on the model of gravitational collapse. The evaporation 
of black holes fl], on which their experimental search 
is based, has not been completely studied either Thus, 
PBH detection will provide valuable information about 
the early Universe and can be a unique test of the general 
theory of relativity, cosmology, and quantum gravity |2 |. 
Direct searches for the bursts of gamma rays from the 
evaporations of PBHs have been carried out in several 
works during the last 15 years |3], H, 0, @, Q, |8|. 

The search for high energy gamma ray bursts on a 
shower array reduces to the search of space and time 
correlations (clusters) of registered extensive air showers 
(EAS). If time intervals used for the search are rather 
small, this search is performed in a horizontal reference 
frame. For each shower i having absolute registration 



time ti and arrival angles {9, (t))i the cluster of events 
i,i + l,i + 2,...,i + n — 1 is sought for, using a 
condition that arrival directions should differ less than 
ar from the weighted mean direction. Thus, each cluster 
is characterized by multiplicity n, duration At, absolute 
time T, and arrival direction {9, (p). 

During such a search, experimentally obtained depen- 
dencies (e.g., cluster registration frequencies for each n) 
are compared to the ones expected from the background 
of accidental coincidences. If measured frequencies of 
cluster registration can be explained with distributions 
expected for accidental coincidences, one can obtain the 
constraints on a gamma ray burst frequency producing 
clusters with a given multiplicity (and thus having par- 
ticular energy flux). 

Evaporating PBHs which have almost reached their 
final evaporation state are a possible source of high 
energy gamma ray bursts. Since the calculated temporal 
and energy characteristics of such bursts depend on the 
theoretical evaporation model |9|, the technique of an 
experimental search and the constraints imposed on the 
PBH number density in the local Universe are model 
dependent. PBHs can be sought for in experiments on 
arrays designed to detect EASs from cosmic rays with 
effective primary gamma-ray energies of 10 TeV or 
higher only within the assumption that the evaporation 
model without a chromosphere |10| is correct (in such 
a model, evaporated particles do not interact with each 
other). The duration of the high-energy GRBs predicted 
by chromospheric evaporation models is too short, much 
shorter than the dead time of EAS arrays, and other 
methods have to be applied for the PBHs search in this 
case 1 11 1. 

It should be noted that the duration of the high energy 
GRBs is fairly short in the evaporation model without a 
chromosphere as well. Therefore, the effect of the array 
dead time on the burst detection probability should be 
taken into account when interpreting the experimental 
data from EAS arrays with a high threshold energy of 
the primary photons. Such analysis have been performed 
in our previous work JS). 

In the experiments searching for gamma ray bursts 
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Fig. 1 : Probabilities of EAS detection by the Andyrchy 
array versus primary gamma-ray photon energy for 
various zenith angles. 



produced by PBHs, the constraints on the PBH number 
density are usually obtained from a condition of absence 
of clusters with a given multiplicity and duration. Thus, 
an one-dimensional task is solved - each cluster (not 
depending on multiplicity) corresponds to the point with 
coordinate Ai on a time axis, and cluster duration 
distribution is analyzed then. However, for the cluster 
with n showers, the number of independent parameters 
(time intervals) is n — 1 and such a cluster can be 
represented as a point in space with n — \ dimensions. 
One can expect that in those cases when time character- 
istics of PBH produced clusters are different from ones 
arising due to accidental coincidences, the increase the 
number of parameters can help to separate PBH events 
from background ones. In this case, due to background 
reduction, the constraint on the PBH number density 
obtained on a given array can be improved. 

This paper is devoted to exploring this idea and 
improving the constraint on PBH number density ob- 
tained previously on Andyrchy array lH). The detailed 
description of this experimental facility can be found in 

inij. 

II. The experiment 

The detection probabilities P{E^, 9) of the EASs gen- 
erated by primary gamma-rays with energy falling 
on Andyrchy array at zenith angle 9 were determined 
by simulating of electromagnetic cascades in the atmo- 
sphere and the detector (see Fig.[I]l |13|. Total number of 
gamma rays which can be detected by the array (integral 
burst profile) is given by 

N^{9,ti) ^ J dE^P{E^,e)dN^/dE^. (1) 



It depends on the time ti left until the end of PBH 
evaporation. 



Fig. 2: Curves 2,3,4: probabilities P(At) for the clusters 
with n = 2, 3, 4 to have the duration in the time interval 
between and At, for zenith angle 9 = 30°. Curve 1 
shows, for comparison, the probability P{ti) for photons 
to be evaporated during time interval between and ti. 

In our previous work fSl we have defined the burst 
duration tb for the given array as time until the end of 
PBH evaporation during which 99% of photons which 
can be detected by this array are evaporated. This 
interval was used as the interval for search of clusters 
from PBHs, i.e., Ts — h (and shower clusters with 
At < tb were sought for). 

Because the exact moment of PBH evaporation is 
unknown, the obtained experimental distributions for 
cluster lengths should be compared not with the burst 
profile, but with the calculated distribution of clusters 
with given multiplicity. Fig. |2] shows, for 9 = 30°, 
the calculated integral At distributions for clusters with 
n — 2,3,4 (curves 2,3 and 4). In the calculation, the 
dead time of the array had been taken into account. In 
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Fig. 3: The dependencies of At9g(n) on zenith angle for 
cluster multiplicity n = 2, 3, 4, 5. 
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Fig. 4: Probabilities of registering the gamma ray clus- 
ters from the evaporating PBH, for the Andyrchy array, 
as a function of zenith angle. 



the same Fig. |2]the burst profile (i.e., the probability to 
detect a gamma particle from evaporating PBH in the 
interval of time between t — ti and full evaporation) is 
also shown (curve 1). It is seen from this figure that 99%- 
level for clusters is achieved at durations much larger 
than <{,. 

Correspondingly, the time interval for search Tg 
should be taken now as = Atgg{n) which we 
define as the time interval in which 99% of clusters 
with multiplicity n fall. Fig. |3] shows the dependence 
of AtgQ{n) on zenith angle for different n. One can 
note that time intervals At99(n) shown in this figure 
are larger than ti, (e.g., for 6 — 0°, ti, = 40 ms - 
see [8 1). Because of this, maximum multiplicities for 
registered clusters also increase: for the case Tg = h, 
maximum multiplicities were 4,4,3,2 for zenith angles 
9 = 0°, 10°, 20°, 30°; for the case Ts = At9g{n) these 
values are 4,4,4,3 for the same corresponding values of 
zenith angles. 

The probability to detect all showers hitting the array 
in a cluster with multiplicity n was calculated using 
Monte-Carlo method taking into account the dead time 
of the array (only those clusters were taken for which 
time interval between consecutive events was larger than 
a dead time of the array <d ^ 1 ms), see Fig. H) 

The result of the search for shower clusters from 
evaporating PBHs is the set of duration distributions of 
clusters for different n and 9. Fig. |5] shows the experi- 
mentally measured distributions for a full observational 
time « 1100 days, summarized over all zenith angles. 
For the search of clusters, we took = 7.0° - such 
a region contains 90% of events from a point source. 
In the same figure, we show modeled distributions for 
clusters produced by accidental coincidences of regular 
EASs. 

It is seen from Fig. |5] that, generally, experimentally 
registered cluster distributions can be explained as acci- 



Fig. 5: Integral duration distributions for clusters with 
multiplicity n = 2, 3, 4, 5, detected by the array for the 
whole region of zenith angles. Points represent the ex- 
perimental results; lines show the calculated distributions 
expected from the cosmic ray background. 



dental coincidences of EASs produced by cosmic rays. 
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Fig. 6: 3-dimensional distributions for clusters with n = 
4, zenith angle is 9 = 30°. a) Events, expected from the 
evaporating PBH; b) events registered by the Andyrchy 
array. 
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III. Multidimensional analysis 

For clusters with n showers, as we have seen, number 
of independent parameters (time intervals) is n — 1, 
and such a cluster can be represented as a point in 
space with corresponding number of dimensions. As an 
example. Fig. |6] shows 3-dimensional distributions of 
clusters expected from PBHs (Fig. |6l-a) and experimen- 
tally registered clusters (Fig. |6}b) for the case of 71 = 4 
and Q = 30°. It is seen that the measured clusters and 
the ones expected from PBHs occupy different regions 
in 3-dimensional space of time intervals. 

In the subsequent analysis, the probabilities for the 
modeled PBH event to get into a (n — 1) -dimensional 
cell with size 1 ms in each dimension were calculated for 
each range of zenith angles 9 and cluster multiplicities 
n. These probabilities were summarized for each cell 
which contained experimentally registered events. If this 
sum did not exceed lO^'^, we assumed, for PBH number 
density constraint calculation, that number of registered 
PBH events is zero (for particular values of n and range 
of zenith angles). We have obtained 3,3,3, and 2 for 
maximum cluster multiplicities of registered events, for 
61 = 0°, 10°, 20°, 30°, correspondingly. 

IV. Results and conclusions 

The constraints on the number density of evaporating 
PBHs in the local region of Galaxy were obtained using 
technique described in work l?). Let a PBH be located 
at distance r from the array and be seen from it at 
zenith angle 9. The mean number of gamma-ray photons 
detected by the array over the burst duration is then 



eNitb{9))Si9) 



(2) 
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where S{9) is the array area and e = 0.9 is the fraction 
of the events from a point source that fell into an angular 
cell with a size of a^. The number of bursts detected 
over the total observation time T can be represented as 

N = PpbhTI/eff , (3) 

with effective volume of the space surveyed by the array 



dfl / drr^F{ni9),n{9)). 



(4) 



Here, F{n, n) is the detection probability of a cluster of 
n or more EASs with the mean value of h. It can be 
expressed through the Poisson probabilities for having i 
showers (and i — n,n + 1, ...) with average value of n 
and probabilities P{i{9),td) [see Fig. HJ to detect all i 
showers having a dead time for one event t^: 



F{n,n) = Y,P{i{9),U) 



(5) 



For the calculation of the effective volume, we take 
n{9) = n'{9) + 1 (this means that the distributions of 
the detected clusters in multiplicity can be explained 
by the background of accidental coincidences; ti'{9) is 



the maximal multiplicity of the detected cluster obtained 
using multidimensional analysis). 

Numerically, we obtain Vcs — 2.8 x 10^^ pc^. If 
the evaporating PBHs are distributed uniformly in the 
local region of the Galaxy, then the upper limit 
on the number density of evaporating PBHs at the 99% 
confidence level can be calculated from the formula 
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As a final result, we obtain pii,„ = 5.4 x 10* pc^^yr^^. 
This is slightly better than the result obtained in our 
previous work |8| with the same set of data but without 
using multidimensional analysis {pum = 8.2 x 10* 
pc~^yr~^). 

The detailed comparison of the results of different 
experiments searching for evaporating PBHs can be 
found in ||8]. This work slightly improves the PBH 
number density constraint of [SJ, but it is still signif- 
icantly weaker than the best (to date) limit obtained 
in the experiment on the Whipple Cherenkov telescope 
[7J (piim = 1-08 x 10^ pc~3yr"i). However, it should 
be noted that the effective gamma-ray energy in our 
experiment is two orders of magnitude higher than that 
of the Whipple telescope. Thus, our upper limit pertains 
not to black holes in general, but to black holes with 
certain properties (those emitting 100-TeV gamma-rays 
at the end of their evaporation during bursts lasting ^10 
ms). 

Acknowledgements. This work was supported by the 
the Russian Foundation for Basic Research (Grants No. 
06-02-16135, 08-07-90400 and 09-02-90900). This work 
was also supported in part by the "Neutrino Physics 
and Astrophysics" Program for Basic Research of the 
Presidium of the Russian Academy of Sciences and 
by "State Program for Support of Leading Scientific 
Schools" (Project No. NSh-321.2008.2). 

References 

[1] S. W. Hawking, Commun. Math. Phys. 43, 199 (1975) [Erratum- 

ibid. 46, 206 (1976)]. 
[2] B. J. Carr, L ect. Notes Phys. 631, 301 (2003) 

|arXiv:astro- ph/03 l0838 |. 
[3] D. E. Alexandreas et al, Phys. Rev. Lett. 71, 2524 (1993). 
[4] M.Amenomori et al. , Proc. 24th International Cosmic Ray 

Conference, Rome (Italy), v.2, 112 (1995). 
[5] B.Funk et al. , Proc. 24th International Cosmic Ray Conference, 

Rome (Italy), v.2, 104 (1995). 
[6] V. Connaughton et al, Astropart. Phys. 8, 179 (1998). 
[7] E. T. Linton et al., JCAP 0601, 013 (2006). 
[8] V. B. Petkov et al, Astron. Lett. 34, 509 (2008) [Pisma Astron. 

Zh. 34, 563 (2008)] rarXiv:0808.3093 [astro-ph]]. 
[9] E. Bugaev, P. Klimai and V. Petkov, Proc. 30th Interna- 
tional Cosmic Ray Conference, Mexico , v.3, 1123 (2007). 

arXiv:0706.3778 [astro-ph]. 
[10] J. H. MacGibbon and B. R. Webber, Phys. Rev. D 41, 3052 

(1990). 

[11] V. B. Petkov, E. V. Bugaev, P A. Klimai and D. V. Smimov, 

JETP Lett. 87, 1 (2008) [Pisma Zh. Eksp. Teor. Fiz. 87, 3 (2008)] 

|ai-Xiv:0803.2313 [astro-ph]]. 
[12] V. B. Petkov et al, Instrum. Exp. Tech. 49, 785 (2006) [Prib. 

Tekh. Eksp. 49, 50 (2006)]. 
[13] D. V. Smimov et al. , Bull. Rus. Acad. Sci.: Physics, 69 N3: 413 

(2005). 



